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Abstract  Genetic  viability  of  threatened  and  endangered 
species  is  of  increasing  concern  with  habitat  loss  and 
fragmentation.  Valuable  assessments  of  the  genetic  status 
of  endangered  species  are  difficult  in  most  cases,  where 
only  single  sample  estimates  are  available.  Using  historical 
and  contemporary  samples,  we  assessed  the  impact  of  both 
historical  and  recent  demographic  changes  on  population 
genetics  of  the  endangered  golden-cheeked  warbler, 
(. Dendroica  chrysoparia).  Our  study  documents  a  steep 
decline  in  genetic  diversity  in  an  endangered  species  over  a 
100-year  period,  along  with  concurrent  increase  in  genetic 
differentiation,  and  low  contemporary  effective  sizes  for  all 
the  populations  we  evaluated.  While  adding  to  the  growing 
body  of  literature  that  describes  the  genetic  impacts  of 
habitat  fragmentation,  our  study  may  also  serve  as  an 
informative  guide  to  future  management  of  endangered 
species.  Our  study  underlines  the  importance  of  long  term 
population  genetic  monitoring  in  understanding  the  full 
extent  of  genetic  changes  in  endangered  species. 
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Habitat  loss  and  fragmentation  are  among  the  greatest 
threats  to  biodiversity  in  the  21st  century  (Young  and 
Clarke  2000).  Species  that  are  habitat  specialists  and 
require  large  home  ranges  are  particularly  sensitive  to 
habitat  loss  and  fragmentation  and  may  experience 
increased  genetic  drift  and  reduced  gene  flow.  While  a 
number  of  studies  have  examined  the  role  of  fragmentation 
on  genetic  structure  (Young  et  al.  1996;  Martinez-Cruz 
et  al.  2007;  Howeth  et  al.  2008;  Sato  and  Harada  2008; 
Yamamoto  et  al.  2004),  these  assessments  are  often 
restricted  to  contemporary  patterns,  when  examination  of 
temporal  changes  would  provide  considerably  more  infor¬ 
mation.  Despite  the  burgeoning  evidence  for  the  implica¬ 
tions  of  habitat  fragmentation  on  contemporary  populations 
(Rodriguez -Munoz  et  al.  2007;  Segelbacher  et  al.  2008),  it 
is  difficult  to  tease  out  the  magnitude  or  trajectory  of 
demographic  and  genetic  change  accompanying  fragmen¬ 
tation,  so  as  to  inform  conservation  management.  This  is 
usually  compounded  by  the  unavailability  of  pre-frag¬ 
mentation  or  ancestral  population  samples  as  reference. 
Although  typing  specimens  from  museum  archives  is  being 
embraced  for  various  evolutionary  and  population  genetic 
studies,  its  application  has  been  relatively  limited,  owing  to 
incomplete  geographical  representation  or  limited  sam¬ 
pling.  This  is  especially  so  for  direct  assessments  of  the 
effects  of  recent  habitat  fragmentation  on  contemporary 
population  structure  and  genetic  diversity.  However,  when 
applied,  this  approach  has  yielded  valuable  demographic  or 
conservation  relevant  information  in  other  species  (Nichols 
et  al.  2001). 
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Identifying  temporal  changes  in  genetic  variation  may  be 
especially  valuable  in  the  case  of  endangered  species,  as 
such  data  might  indicate  whether  the  effective  population 
size  has  been  small  enough  to  result  in  reduced  genetic 
variation  (and  thus  adaptive  potential)  and  if  populations 
have  become  recently  isolated  and  thus  subject  to  local 
extirpation.  Temporal  changes  in  genetic  variance  are  also 
valuable  for  understanding  the  effective  size  ( Ne )  in  relation 
to  the  census  size  ( Nc )  because  of  the  disparity  between  the 
number  of  individuals  in  a  population  and  actual  genetic 
contributors  to  each  generation  (Nei  and  Tajima  1981; 
Crow  and  Denniston  1988;  Frankham  1995a).  Another 
advantage  of  a  temporal  analysis  is  in  facilitating  quantifi¬ 
cation  of  genetic  change  in  populations,  partly  satisfying  the 
proximate  data  demands  in  a  conservation  context. 

We  present  our  results  from  temporal  analyses  of 
genetic  data  from  the  endangered  golden-cheeked  warbler 
( Dendroica  chrysoparia)  from  samples  spanning  the  years 
1890-2005.  In  this  study  we  were  able  to  utilize  archived 
museum  samples  as  reference  populations  to  understand 
the  impacts  of  habitat  fragmentation  on  contemporary 
populations  of  golden-cheeked  warblers.  This  migratory 
songbird  has  a  breeding  range  restricted  to  central  Texas, 
USA  (Fig.  1).  Estimates  of  population  sizes  based  on 
banding  data  and  surveys  were  18,500  and  13,800  in  1962 
and  1990,  respectively,  with  a  35%  decline  in  breeding 
habitat  (USFWS  1992).  Preferred  nesting  habitat  for  this 
species  are  climax  forests  of  Ashe  juniper  (Juniperus  ashei) 
and  oak  ( Quercus  spp.)  that  develop  on  soils  derived  from 
limestone  substrates  (Kroll  1980;  USFWS  1992).  Similar 
to  several  other  neotropical  songbirds,  golden-cheeked 
warblers  exhibit  site  fidelity,  territoriality  and  natal 
philopatry — all  of  which  may  have  implications  for  popu¬ 
lation  structure.  Fragmentation  of  available  habitat  is 
known  to  have  resulted  from  clearing  of  woodlands  for 
livestock  production  and  agriculture.  Populations  are  also 
depressed  due  to  brood  parasitism  by  brown-headed  cow- 
birds,  Molothrus  ater  (Sperry  et  al.  2008). 

In  an  earlier  assessment  of  genetic  structure  in  D.  chrys¬ 
oparia  we  found  that  few  populations  were  significantly 
differentiated,  with  little  evidence  of  genetic  bottlenecks, 
and  genetic  diversity  comparable  to  more  abundant  warbler 
species  (Lindsay  et  al.  2008).  As  with  similar  studies,  that 
investigation  was  based  on  contemporary  samples  alone. 
However,  extrapolating  from  data  based  on  a  single  time 
point  limits  the  resolution  of  the  full  extent  of  recent  and 
ongoing  genetic  changes.  More  complete  assessments  of 
genetic  changes  require  assessing  the  temporal  component  of 
genetic  change.  Here  we  used  samples  taken  from  2  to  3  time 
points  spanning  across  a  century,  for  three  sites  in  central 
Texas  (Fig.  1)  to  examine  temporal  change  in  genetic  char¬ 
acteristics  and  estimate  the  effective  population  size  ( Ne ). 


Materials  and  methods 

Locations  and  times  of  samples  for  temporal  analyses 

Our  study  focused  on  temporal  samples  collected  from 
multiple  locations  between  the  years  1890  and  2008.  This 
sampling  window  provided  the  greatest  extent  of  temporal 
separation  between  samples  and  also  brackets  the  period 
during  which  extensive  habitat  loss  and  fragmentation  is 
documented  to  have  taken  place.  Due  to  limited  archived 
material  available  in  museum  collections,  sampled  loca¬ 
tions  correspond  to  only  three  (Fig.  1 )  of  seven  populations 
studied  in  Lindsay  et  al.  (2008).  We  categorized  compari¬ 
sons  of  genetic  diversity  and  effective  population  size 
within  and  between  time  points  into  two  classes — historical 
sample  comparisons,  and  historical-contemporary  com¬ 
parisons  (Table  1).  We  refer  to  pre-decline  samples  from 
1890  to  1915  as  the  historical  samples,  while  samples  from 
2005  are  referred  to  as  the  contemporary  samples.  Histor¬ 
ical  sample  locations  correspond  to  the  same  county  from 
which  contemporary  samples  were  collected  except  for 
Bexar,  which  comprised  samples  from  Bexar-Comal- 
Kendall  (<50  km  radius).  The  contemporary  Bexar  and 
Travis  samples  correspond  to  the  CAB  and  BAC  samples 
from  Lindsay  et  al.  (2008),  respectively.  Samples  labeled 
1900  for  Kerr  consist  of  samples  collected  1895-1902, 
whereas  samples  labeled  1890  for  Bexar  consist  of  samples 
collected  1885-1895. 

Archived  samples  and  processing 

We  followed  the  guidelines  of  Bonin  et  al.  (2004)  in  pro¬ 
cessing  of  historical  and  contemporary  samples  to  mini¬ 
mize  genotyping  errors.  Historical  samples  were  obtained 
from  four  national  ornithological  collections,  namely  the 
National  Museum  of  Natural  History  (NMNH),  Washing¬ 
ton,  DC,  USA,  the  Museum  of  Comparative  Zoology 
(MCZ),  Cambridge,  MA,  the  Field  Museum  of  Natural 
History  (FMNH),  Chicago,  IL  and  the  American  Museum 
of  Natural  History  (AMNH),  New  York,  NY  (sample  IDs 
are  provided  in  Appendix  2).  A  total  of  73  specimens  were 
sampled  and  from  each  archived  specimen,  a  small  piece 
(1  mm3)  of  skin  was  extracted  from  the  toe-pads  (Ellegren 
1991)  with  a  sterile  razor  blade  and  stored  in  a  sterile 
microcentrifuge  tube.  DNA  extraction  was  performed 
using  a  QiaAmp  Microextraction  kit  (Qiagen,  Valencia, 
CA),  with  initial  overnight  digest,  followed  by  spin-column 
DNA  isolation  and  final  elution  in  40  pi  low-TE  buffer. 
Samples  were  stored  in  low-TE  buffer  at  — 20°C  to  prevent 
degradation  of  potentially  fragile,  low-yield  DNA  extracts 
until  genotyping.  All  archived  sample  processing  was 
performed  in  a  separate  lab  under  a  sterile  hood. 
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Fig.  1  Map  showing  the 
distribution  of  breeding  habitat 
for  the  D.  chrysoparia  in  central 
Texas.  Though  the  shaded 
region  covers  the  extent  of  the 
current  breeding  range, 
available  breeding  habitat  is 
actually  relatively  rare  and  is 
highly  fragmented.  Filled  ovals 
represent  sites  where  both 
historical  and  contemporary 
samples  have  been  procured 


Table  1  Sample  locations,  time  points  and  sample  sizes  used  for 
temporal  analysis  of  changes  in  genetic  diversity,  population  differ¬ 
entiation  and  estimation  of  the  effective  population  size  Ne 


Sampling  site 

Sample  year 

Sample  size  ( N ) 

Bexar 

1890 

11 

2005 

17 

Kerr 

1900 

9 

1915 

33 

2005 

31 

Travis 

1913 

16 

2005 

17 

Establishing  reliability  and  estimating  error  rates 

We  first  established  a  framework  to  quantify  the  reliability 
of  our  genotype  dataset  using  the  program  RELIOTYPE 
(Miller  et  al.  2002).  This  program  uses  a  maximum-like¬ 
lihood  approach  for  reducing  genotyping  errors  and  pro¬ 
vides  recommendations  for  replications  at  each  locus  to 
reach  a  >95%  reliability  level  (the  minimum  acceptable 
level).  Once  individual  loci-sample  combinations  reach  or 


exceed  the  reliability  threshold,  the  data  can  be  considered 
to  have  a  negligible  rate  of  error.  We  initially  determined 
consistent  amplification  if  at  least  three  out  of  four  attempts 
resulted  in  the  expected  Polymerase  Chain  Reaction  (PCR) 
product.  We  used  MICROCHECKER  (Van  Oosterhout 
et  al.  2004)  to  check  our  data  for  genotyping  errors,  large- 
allele  dropout  and  null  alleles.  For  loci  that  passed  these 
initial  criteria,  we  analyzed  the  data  in  RELIOTYPE  to 
determine  if  and  how  many  replicates  were  necessary  to 
meet  a  minimum  reliability  level  of  95%  (for  each  locus). 
Amplification  of  loci  that  did  not  meet  the  reliability 
threshold  were  replicated  for  the  recommended  number  of 
times  and  analyzed  again  to  confirm  that  they  exceeded  the 
95%  threshold.  In  the  final  dataset,  we  retained  only  loci 
with  average  reliability  scores  equaling  or  greater  than 
98%.  We  also  used  the  program  GIMLET  (Valiere  2002)  to 
estimate  the  error  rate  arising  due  to  allelic  dropout  across 
all  replicates. 

Each  sample  was  genotyped  following  amplification  by 
PCR  at  each  of  nine  microsatellite  loci  previously  char¬ 
acterized  for  this  species  (Lindsay  et  al.  2008),  followed  by 
fragment  analysis  on  an  ABI  3130  Genetic  Analyzer.  PCR 
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optimization  and  genotyping  analyses  were  performed  as 
reported  in  Lindsay  et  al.  (2008).  PCR  on  DNA  from 
archived  samples  was  optimized  by  changes  to  MgCl2 
concentrations  or  with  thermoprofile  variations. 

Statistical  analysis 

We  collected  allelic  data  following  fragment  analyses 
using  GENESCAN  software  (version  3.1,  Applied  Bio¬ 
systems,  Foster  City,  CA).  We  performed  exact  tests  for 
deviations  from  Hardy- Weinberg  expectations  on  allele 
frequencies  using  the  online  tool  GENEPOP  ON  THE 
WEB  (Raymond  and  Rousset  1995).  The  program  GE- 
NETIX  (Belkhir  et  al.  2004)  was  used  to  estimate  measures 
of  unbiased  expected  heterozygosity  (HEXp),  and  the  pro¬ 
gram  FSTAT  (Goudet  2002)  was  used  to  estimate  allelic 
richness  ( AR ).  AR  controls  for  biases  in  allelic  diversity 
arising  from  unequal  sample  sizes  (Leberg  2002).  Global 
diversity  indices  and  other  sample-by-locus  statistics  are 
provided  in  Appendix  1. 

To  quantify  the  changes  in  genetic  diversity  between  the 
historical  and  contemporary  samples,  we  used  a  paired 
/-test  (PROC  TTEST,  Statistical  Analysis  Software,  SAS 
Institute,  Cary,  NC),  by  locus,  to  compare  HEXp  and  AR 
between  time  points. 

Historical  population  genetic  structure  was  characterized 
by  estimating  pairwise  FST  among  the  historical  samples 
from  Kerr,  Bexar,  and  Travis  populations  in  GENETIX 
(Belkhir  et  al.  2004).  The  significance  of  E-statistic  esti¬ 
mates  was  assessed  using  a  resampling  approach  (3,000 
permutations)  to  account  for  differences  in  population 
sample  sizes.  Identical  analyses  were  performed  on  the 
corresponding  contemporary  populations.  To  verify  that 
obtained  F st  estimates  are  not  an  artifact  of  high  intra¬ 
population  diversity,  we  also  compared  differentiation 
using  two  other  estimator,  namely  G’st  and  DEST  (Hedrick 
2005;  Jost  2008).  We  used  a  paired  (by  locus)  /-test  to 
determine  if  the  degree  of  pairwise-differentiation  was 
different  between  the  historical  and  contemporary  periods. 
Additionally,  elevated  linkage  disequilibrium  has  been 
associated  with  genetic  fragmentation  (Zartman  et  al.  2006) 
and  we  tested  whether  the  differentiation  we  observed  also 
caused  elevated  linkage  disequilibrium. 

Estimation  of  Ne 

We  used  the  pseudo-maximum  likelihood  estimator  (Wang 
2001;  Wang  and  Whitlock  2003)  as  implemented  in  the 
program  MLNE  to  estimate  Ne  (hereafter  referred  to  as 
ML/Ve).  MLNE  provides  estimates  of  the  variance  Ne 
across  time- spans  separating  two  or  more  sampling  inter¬ 
vals  (Wang  and  Whitlock  2003).  We  estimated  Ne  jointly 
with  migration  rate  (m),  where  it  is  assumed  that  migrants 


originate  from  an  infinitely  large  population.  Estimation  of 
m  requires  that  a  source  population  be  specified  from  which 
gene  flow  may  occur.  For  each  of  the  three  populations  for 
which  we  estimated  ML/Ve,  we  specified  a  source  popula¬ 
tion  comprised  of  samples  combined  from  the  remaining 
two  populations.  We  used  this  estimator  as  it  is  the  only 
available  temporal  estimator  that  considers  populations 
open  to  migration,  and  as  we  were  specifically  interested  in 
the  temporal  estimate  of  Ne.  For  comparison  with  the 
M  LNe  estimates,  we  also  estimated  the  moment  Ne  (Mt Ne) 
based  on  equations  15  and  18  from  Nei  and  Tajima  (Nei 
and  Tajima  1981),  which  assumes  closed  populations. 

Estimation  of  Ne  using  either  the  moment  or  likelihood 
approach  requires  an  assumption  of  discrete  generations. 
As  D.  chrysoparici  do  not  conform  to  the  assumption  of 
discrete  generations,  we  used  two  different  generation 
lengths  (T  =  1  and  2  years)  to  evaluate  the  effects  that 
T  may  have  on  estimates  of  Ne.  The  majority  of  adult 
females  in  songbirds  breed  in  their  first  year,  whereas 
reports  indicate  that  age  of  males  at  first  breeding  may  be 
either  one  or  two  year  (USFWS  1992).  Considering  this, 
either  generation  length  may  be  appropriate,  and  we 
applied  both  in  separate  analyses  in  addition  to  the  average 
of  the  two  (1.5  year/generation).  M LNe  estimates  obtained 
from  T  =  1  and  T  —  2  were  different,  but  were  qualita¬ 
tively  similar,  with  results  for  T  =  1.5  being  intermediate 
of  the  two.  Therefore,  only  results  from  T  —  1  are  reported 
in  the  main  text  (Table  3).  Estimates  for  T  —  2  are  pre¬ 
sented  in  Table  SI. 

Current  census  population  sizes  ( Nc )  for  the  three  sites 
are  approximately  90  (Frels  Jr  2006),  485  (Source:  Camp 
Bullis,  TX)  and  1,000  (USFWS  2004)  pairs,  for  Kerr, 
Bexar,  and  Travis,  respectively  (Lindsay  et  al.  2008). 
Census  methodologies  vary  between  sites  and  are  subject 
to  minor  yearly  variations,  hence  population  sizes  are  not 
to  be  considered  to  be  exact  and  are  reported  as  3-year 
averages. 

Results 

Approximately  78%  of  PCR  amplifications  on  archived 
samples  amplified  on  the  first  attempt,  compared  to  98% 
for  the  contemporary  samples.  Consensus  genotypes  from 
historical  samples  were  obtained  following  multiple  repli¬ 
cations  (up  to  six  replicates  for  different  loci)  based  on  the 
recommendations  from  the  program  RELIOTYPE.  Only 
eight  of  the  nine  loci  analyzed  were  amplified  consistently 
and  met  the  98%  reliability  threshold  that  we  had  estab¬ 
lished.  The  eight  polymorphic  microsatellite  loci  (Lindsay 
et  al.  2008)  were  scored  in  both  historical  and  contempo¬ 
rary  samples,  with  no  evidence  of  null  alleles.  Error  rates  in 
historical  sample  genotypes  were  estimated  to  be  1.9% 
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across  loci.  A  total  of  95  different  alleles  were  detected  in 
the  historical  samples,  of  which  92  were  found  in  con¬ 
temporary  samples.  We  did  not  observe  significant  devia¬ 
tions  from  Hardy- Weinberg  equilibria  across  loci. 

Based  on  our  analyses  of  temporal  changes  in  genetic 
variation  over  a  100-year  interval,  we  found  that  popula¬ 
tions  of  the  D.  chrysoparia  show  reduced  genetic  diversity, 
increased  genetic  fragmentation  and  small  effective  popu¬ 
lation  sizes  compared  to  historical  estimates.  There  was  no 
difference  between  1900  and  1915  estimates  of  allelic 
richness  (AR)  and  heterozygosity  ( HEXp )  (Fig.  2). 

Both  Ar  and  HEXP  were  significantly  lower  in  contem¬ 
porary  samples  compared  to  historical  samples  from  the 
same  sites,  with  average  declines  of  20  and  13%,  respec¬ 
tively  (Fig.  2).  Though  the  two  estimates  showed  similar 
declines,  there  was  a  greater  proportional  change  in  AR. 

No  significant  genetic  differentiation  was  detected 
among  the  three  historical  populations  in  our  sample  as 
measured  by  f'si  (Table  2),  and  we  recovered  a  similar, 
albeit  slightly  different,  pattern  of  differentiation  among 
the  three  contemporary  populations  as  reported  in  Lindsay 
et  al.  (2008).  The  most  important  result,  however,  was  that 
measures  of  pairwise  differentiation  had  increased  signifi¬ 
cantly  when  comparing  historic  to  contemporary  levels 
of  differentiation  for  the  same  population  pairs.  This  result 
was  retained  with  other  estimators  of  differentiation 
( G’ ST '>  Dst),  which  control  for  high  genetic  diversity 
within  populations  (see  Appendix  1).  We  also  found  that 
~  56%  of  pairwise  locus  tests  were  significant  for  linkage 
disequilibrium. 

Across  our  comparisons,  MtNe  almost  always  produced 
estimates  larger  than  ML/Ve  (Table  3);  Mt/V,,  estimates  of 
historical  Ne  for  Kerr  (1900-1915)  were  very  large  (co)  in 
all  cases  (Table  3),  whereas  the  ML/Ve  =  144  (95%  Cl 
39-6242).  The  estimated  rate  of  migration  was  also  high  in 
proportion  to  Ne  (Table  3). 

The  historical-contemporary  estimates  of  Ne  for  Bexar 
(1890-2005),  Travis  (1913-2005),  and  Kerr  (1915-2005) 
were  49  (95%  Cl  24-212),  273  (95%  Cl  46-621),  and  66 
(95%  Cl  43-419),  respectively  (Table  3).  As  our  estimates 
of  Fst  indicate  that  rates  of  gene  flow  may  not  have  been 
constant  over  the  period  separating  the  samples,  we  are  not 
reporting  estimates  of  m. 

Discussion 

Understanding  the  genetic  impacts  of  fragmentation  is 
clearly  of  interest  for  conservation  of  species,  and  the 
increasing  attempts  to  quantify  these  impacts  in  various 
taxa  attest  to  the  exigency  of  the  threats  faced  by  biodi¬ 
versity  (Alo  and  Turner  2005;  Veit  et  al.  2005;  Martinez- 
Cruz  et  al.  2007;  Born  et  al.  2008;  Jordan  and  Snell  2008; 


Time 

Fig.  2  Estimates  of  mean  genetic  diversity  (with  SEs)  for  historical 
and  contemporary  populations  of  D.  chrysoparia.  The  hypotheses  of 
no  difference  in  estimates  of  allelic  richness  (Afi)  and  expected 
heterozygosity  (HEXp)  were  evaluated  with  a  T-test  (a  =  0.05). 
Symbols  above  each  bar  indicate  significant  pairwise  differences 
(*P  <  0.05;  **P  <  0.01;  ***P  <  0.001) 


Kramer  et  al.  2008;  Angelone  and  Holderegger  2009). 
Despite  these  efforts,  the  evidence  for  the  genetic  impacts 
of  fragmentation  is  not  unequivocal.  Our  study  addresses 
some  of  the  lingering  questions  by  analyzing  genetic 
changes  from  samples  straddling  a  period  of  well  docu¬ 
mented  habitat  loss  and  fragmentation  in  central  Texas. 
Our  results  follow  similar  reports  of  reduced  genetic 
diversity,  increased  spatial  structuring  or  low  effective 
population  sizes  in  various  taxa,  resulting  from  habitat 
fragmentation  (Ellis  et  al.  2006;  Chevolot  et  al.  2008;  Dixo 
et  al.  2009).  However,  our  study  is  among  the  few  that  were 
able  to  report  a  decline  by  comparison  with  pre-fragmen¬ 
tation  samples.  Our  results  add  to  the  important  body  of 
literature  on  the  effects  of  fragmentation  on  conservation 
genetics  and  are  potentially  informative  for  future  man¬ 
agement  of  the  golden-cheeked  warbler. 

Changes  in  diversity 

Large  losses  of  AR  are  characteristic  of  bottlenecks  (Leberg 
1992)  and  are  of  concern  because  random  loss  of  alleles 
might  reduce  adaptive  potential.  Although  the  decline  in 
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Table  2  Pairwise  differentiation  ( Fst )  between  pairs  of  sites  that  were  sampled  for  both  historical  and  contemporary  periods 


Sample  pair 

Historical  samples 

Contemporary  samples 

Historical  v/s  contemporary  comparison 

Fst 

P 

Fst 

P 

T 

P 

Bexar-Kerr 

0.00053 

0.485 

0.031 

<0.01 

16.46 

<0.001 

Travis-Bexar 

0.00093 

0.436 

0.051 

<0.01 

12.85 

<0.001 

Travis-Kerr 

0.000095 

0.455 

0.009 

0.05 

19.40 

<0.001 

Within  each  time  period,  we  tested  the  hypothesis  that  pairwise  FST  estimates  would  not  be  significantly  different  from  zero,  the  expectation  for 
panmictic  populations.  For  temporal  comparisons,  we  used  paired  T-tests  to  test  the  hypothesis  that  pairwise  Fst  estimates  among  populations  did 
not  differ  between  historical  and  the  contemporary  periods 


Table  3  Estimation  of  effective  population  sizes  (Ne)  using  temporal  sampling  approaches  and  a  generation  length  of  1  year 


Likelihood  estimator 


Site 

Time 

No.  gen 

MfiVe 

MLW, 

95%  Cl 

m 

95%  Cl 

Nc 

Historical  sample  comparisons 

Kerr 

1900-1915 

15 

oo 

144 

39-6242 

0.362 

0.016-0.899 

- 

Historical— 

-contemporary  comparisons 

Bexar 

1890-2005 

115 

98 

49 

24-212 

0.084 

0.019-0.187 

970 

Kerr 

1915-2005 

90 

326 

66 

35-725 

0.053 

0.012-0.365 

184 

Travis 

1913-2005 

92 

378 

273 

46-621 

0.047 

0.014-0.323 

2000 

Values  (with  95%  confidence  intervals)  include  Nei  and  Tajima’s  moment  estimate  (MtA(,),  Wang’s  likelihood  estimate  (MLlVe),  migration  (m), 
contemporary  census  size  (Nc).  Nc  is  not  available  for  historical  Kerr  populations 


heterozygosity  was  not  to  the  same  extent  as  allelic 
richness,  these  two  indices  change  at  different  rates  in 
response  to  demographic  changes,  suggesting  that  popu¬ 
lations  sizes  may  have  been  low  for  a  sufficient  period  of 
time  for  inbreeding  to  accumulate  (Nei  et  al.  1975; 
Frankham  1995c)  in  D.  chrysoparia  populations.  Losses 
of  A R  are  expected  to  occur  rapidly  following  bottlenecks 
compared  to  changes  in  HEXp  (Young  et  al.  1996)  Recent 
population  estimates  of  D.  chrysoparia  populations  sug¬ 
gest  that  declines  may  have  slowed,  as  a  result  of  which 
declines  in  AR  may  have  ceased,  although  it  is  possible 
that  losses  of  heterozygosity  are  ongoing,  given  the 
slower  rate  at  which  heterozygosity  responds  to  popula¬ 
tion  size  changes  (Nei  et  al.  1975).  Similar  declines  in 
genetic  diversity  have  been  documented  in  other  species 
(Ellis  et  al.  2006;  Lage  and  Kornfield  2006),  but  establish¬ 
ing  a  link  with  fragmentation  has  not  always  been  straight¬ 
forward.  Perhaps  the  most  striking  finding  of  our  study 
is  its  occurrence  in  combination  with  increased  genetic 
differentiation. 

Increased  differentiation 

We  obtained  slightly  different  results  compared  to  Lind¬ 
say  et  al.  (2008)  partly  as  a  result  of  increased  sample 
size  for  the  Kerr  population  (31  instead  of  16)  and  using 


one  fewer  loci  than  in  the  previous  study.  Additionally, 
all  samples  were  reanalyzed  and  scored  on  a  different 
genotyping  platform  compared  to  Lindsay  et  al.  (2008).  It 
is  also  possible  that  the  partitioning  of  variation  would  be 
different  when  using  different  numbers  of  populations  in 
the  analysis  (a  total  seven  populations  were  surveyed 
in  that  case).  Our  results  suggest  reduced  gene  flow 
in  contemporary  populations  compared  to  historical 
times  indicating  that  barriers  to  gene  flow  are  recent 
phenomena.  Additionally,  our  finding  of  elevated  linkage 
disequilibrium,  providing  additional  evidence  for  frag¬ 
mentation.  Historical  accounts  of  habitat  from  central 
Texas  at  turn  of  the  20th  century  indicate  that  preferred 
D.  chrysoparia  habitat  may  have  been  more  contiguous 
than  today  (Schmidly  and  Sansom  2002),  and  rapid 
urbanization  of  central  Texas  since  the  1940s  (including 
225,000  km  of  road  construction)  is  known  to  have 
greatly  reduced  the  extent  and  quality  of  this  songbird’s 
habitat  (USFWS  1992).  In  a  recent  review  on  the  genetics 
impacts  of  landscape  alteration  we  report  that  changes  in 
population  structure  appear  to  be  more  common  in  frag¬ 
mented  populations  than  losses  of  genetic  diversity 
(Leberg  et  al.  2009).  The  finding  of  increased  differenti¬ 
ation  along  with  losses  of  diversity  is  one  of  the  main 
findings  of  our  study,  and  was  achieved  mainly  by  the 
availability  of  temporal  samples. 
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Effective  size  estimates 

We  found  relatively  low  estimates  of  Ne  based  on  two 
different  estimators,  but  they  are  difficult  to  interpret 
without  accurate  knowledge  of  census  size.  The  relation¬ 
ship  between  Ne  and  Nc  is  much  debated  in  literature 
(Frankham  1995b;  Nunney  1995),  but  an  Ne  that  is  sig¬ 
nificantly  lower  than  Nc  indicates  that  a  bottleneck  may 
have  occurred  and  may  negatively  impact  genetic  variation 
over  time;  thus  this  ratio  is  of  obvious  importance  to  spe¬ 
cies  conservation  (Nunney  2000).  We  did  not  estimate  NJ 
Nc  ratios  in  I).  chrysoparia  as  there  is  little  information 
about  historical  census  sizes,  making  such  a  comparison 
impossible.  The  relatively  small  estimates  of  Ne  observed 
in  our  study  may  be  a  consequence  of  various  demographic 
factors  related  to  habitat  loss  and  fragmentation,  and  may 
reflect  poor  survival  in  habitat  edges  (Lindsay  et  al.  2008), 
reduced  nest  success  in  degraded  habitats  (Peak  2007; 
Reidy  et  al.  2009),  high  levels  of  brood  parasitism  and  nest 
predation  (Sperry  et  al.  2008),  or  limitation  of  habitat  on 
wintering  grounds  (Rappole  et  al.  2003).  A  follow-up 
temporal  analysis  of  the  managed  populations  will  be 
required  to  assess  if  and  how  Ne  is  changing  in  response  to 
management  efforts. 

Alternatively,  D.  chrysoparia  populations  may  not  have 
been  very  large  historically.  Although  forested  areas  are 
presumed  to  have  been  more  contiguous  in  historical  times 
than  at  present,  it  is  debatable  whether  suitable  habitat  was 
extensive  (given  the  highly  specialized  habitat  require¬ 
ments).  In  a  scenario  where  D.  chrysoparia  would  have 
historically  existed  at  relatively  small  effective  sizes, 
continuing  fragmentation  might  significantly  impact 
genetic  diversity.  The  observed  steep  increase  in  genetic 
differentiation  indicates  that  fragmentation  becomes  criti¬ 
cal  for  populations  that  are  not  large  to  begin  with,  and 


suffer  further  declines  due  to  habitat  loss.  This  suggests 
that  gene  flow  may  be  much  more  sensitive  to  habitat 
changes  than  to  within  population  demographic  events 
(Leberg  et  al.  2009).  This  apparent  loss  of  connectedness 
among  populations  is  expected  when  habitat  becomes 
fragmented  and  remnant  populations  become  isolated  on 
patches  of  preserved  habitat — the  current  situation  with 
D.  chrysoparia. 

Our  findings  of  a  significant  decline  in  genetic  diversity, 
and  increases  in  genetic  differentiation  are  an  important 
addition  to  the  data  on  fragmentation  effects  on  vertebrate 
species.  Given  the  high  degree  of  habitat  specialization  of 
D.  chrysoparia  and  the  slow  recovery  of  habitat  to  climax 
oak-juniper  communities,  our  results  elicit  additional  con¬ 
cern  for  the  species’  long-term  viability.  Considering  that 
our  earlier  study  on  contemporary  populations  failed  to 
detect  many  of  the  genetic  changes  that  have  occurred  in 
the  last  century,  this  study  emphasizes  the  importance  of 
temporal  sampling  for  characterizing  genetic  trends  in 
endangered  species. 
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Appendix  1 

See  Tables  4,  5,  6  and  7 


Table  4  Historical  population  summary  of  diversity  indices 


Gene  diversity  per  locus  and  population 


Locus  1 

0.945 

0.954 

0.941 

Locus2 

0.898 

0.881 

0.927 

Locus3 

0.852 

0.854 

0.877 

Locus4 

0.846 

0.827 

0.886 

Locus5 

0.873 

0.846 

0.923 

Locus6 

0.867 

0.919 

0.891 

Locus7 

0.866 

0.848 

0.877 

Locus8 

0.865 

0.879 

0.905 

0.8765 

0.876 

0.903 
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Table  4  continued 

Number  of  alleles  sampled 

Locus  1 

19 

18 

12 

22 

Locus2 

11 

10 

12 

12 

Locus3 

9 

9 

7 

9 

Locus4 

9 

9 

8 

10 

Locus5 

10 

8 

10 

10 

Locus6 

11 

11 

8 

11 

Locus7 

12 

10 

9 

14 

Locus8 

13 

10 

10 

16 

11.75 

10.625 

9.5 

Allelic  richness  per  locus  and  population  based  on  min.  sample  size  of:  1 1  diploid  individuals 

Locus  1 

13.084  14.283  12 

13.063 

Locus2 

9.185  9.047  12 

9.359 

Locus3 

7.202  7.875  7 

7.232 

Locus4 

7.491  7.851  8 

7.573 

Locus5 

8.282  7.236  10 

8.148 

Locus6 

8.279  9.885  8 

8.379 

Locus7 

8.154  8.123  9 

8.111 

Locus8 

9.003  8.845  10 

9.196 

8.835  9.1431  9.5 

Table  5  Nei’s  estimation  of  heterozygosity 

LocName 

Ho 

Hs 

Ht 

Dst 

Dst" 

Ht’ 

Gst 

Gst’ 

Gis 

Locus  1 

0.888 

0.947 

0.945 

-0.002 

-0.003 

0.944 

-0.002 

-0.004 

0.062 

Locus2 

0.855 

0.902 

0.902 

0 

0 

0.902 

0 

0 

0.052 

Locus3 

0.831 

0.861 

0.862 

0 

0.001 

0.862 

0 

0.001 

0.035 

Locus4 

0.917 

0.853 

0.853 

0.001 

0.001 

0.853 

0.001 

0.001 

-0.076 

Locus5 

0.868 

0.88 

0.874 

-0.007 

-0.01 

0.87 

-0.007 

-0.011 

0.014 

Locus6 

0.883 

0.892 

0.883 

-0.009 

-0.013 

0.878 

-0.01 

-0.015 

0.009 

Locus7 

0.917 

0.864 

0.854 

-0.009 

-0.014 

0.85 

-0.011 

-0.017 

-0.062 

Locus8 

0.912 

0.882 

0.878 

-0.005 

-0.007 

0.875 

-0.005 

-0.008 

-0.034 

Overall 

0.884 

0.885 

0.881 

-0.004 

-0.006 

0.879 

-0.004 

-0.007 

0.001 

Table  6  Modem  population's  summary  of  diversity 

Gene  diversity  per  locus  and  population 

Loci 

0.93 

0.869 

0.956 

Loc2 

0.759 

0.745 

0.807 

Loc3 

0.721 

0.765 

0.704 

Loc4 

0.532 

0.587 

0.563 

Loc5 

0.829 

0.823 

0.836 
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Table  6  continued 

Gene  diversity  per  locus  and  population 

Loc6 

0.841 

0.752 

0.733 

Loc7 

0.841 

0.868 

0.785 

Loc8 

0.822 

0.843 

0.893 

0.784 

0.7815 

0.7846 

Number  of  alleles  sampled 

Loci 

16 

19 

19 

Loc2 

7 

8 

6 

Loc3 

6 

6 

6 

Loc4 

5 

6 

5 

Loc5 

8 

10 

8 

Loc6 

8 

9 

9 

Loc7 

7 

11 

7 

Loc8 

8 

9 

11 

8.125 

9.75 

8.875 

Allelic  richness  per  locus  and  population  based  on  min.  sample  size 

of:  17  diploid  individuals 

Loci 

15.228 

14.511 

19 

Loc2 

6.789 

6.902 

6 

Loc3 

5.878 

5.461 

6 

Loc4 

4.789 

5.059 

5 

Loc5 

7.789 

8.319 

8 

Loc6 

7.789 

7.72 

9 

Loc7 

6.895 

9.315 

7 

Loc8 

7.877 

8.458 

11 

7.877 

8.2181 

8.875 

Table  7  Diversity  measures 

Locus 

Ho 

Hs 

Ht 

Dst 

Dst’ 

Ht’ 

Gst 

Gst' 

Gis 

Loci 

0.898 

0.918 

0.938 

0.02 

0.03 

0.948 

0.021 

0.031 

0.022 

Loc2 

0.79 

0.77 

0.768 

-0.002 

-0.003 

0.767 

-0.002 

-0.004 

-0.025 

Loc3 

0.813 

0.73 

0.767 

0.037 

0.055 

0.785 

0.048 

0.07 

-0.113 

Loc4 

0.581 

0.561 

0.566 

0.005 

0.008 

0.569 

0.01 

0.014 

-0.035 

Loc5 

0.852 

0.829 

0.827 

-0.002 

-0.003 

0.827 

-0.002 

-0.003 

-0.028 

Loc6 

0.77 

0.775 

0.808 

0.032 

0.048 

0.824 

0.04 

0.059 

0.007 

Loc7 

0.813 

0.831 

0.872 

0.04 

0.06 

0.892 

0.046 

0.068 

0.022 

Loc8 

0.87 

0.853 

0.855 

0.003 

0.004 

0.856 

0.003 

0.005 

-0.02 

Overall 

0.798 

0.783 

0.8 

0.017 

0.025 

0.808 

0.021 

0.031 

-0.02 
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Appendix  2 


Sample  ID 

Museum 

A3  82291 

AMNH 

A3  82292 

AMNH 

A382293 

AMNH 

A3  82294 

AMNH 

A3  82295 

AMNH 

A3  82296 

AMNH 

A3  82297 

AMNH 

A382360 

AMNH 

A382361 

AMNH 

A382362 

AMNH 

A382363 

AMNH 

A382364 

AMNH 

A382365 

AMNH 

A3  82305 

AMNH 

A3  823 10 

AMNH 

A3  823 13 

AMNH 

A3  823 18 

AMNH 

A382323 

AMNH 

A382327 

AMNH 

A382329 

AMNH 

A382335 

AMNH 

A3  82340 

AMNH 

A3  82341 

AMNH 

A3  82346 

AMNH 

A3  82348 

AMNH 

A382351 

AMNH 

A382352 

AMNH 

A507004 

AMNH 

F149870 

FMNH 

F149871 

FMNH 

F149872 

FMNH 

F26172 

FMNH 

F308240 

FMNH 

F308244 

FMNH 

M101045 

MCZ 

M101254 

MCZ 

M204935 

MCZ 

M2 143 64 

MCZ 

M28413 

MCZ 

M326084 

MCZ 

M326086 

MCZ 

M320914 

MCZ 

M320923 

MCZ 

M320940 

MCZ 

M3 20941 

MCZ 

M320943 

MCZ 

M3 20946 

MCZ 

Appendix  continued 


Sample  ID 

Museum 

M320956 

MCZ 

M320957 

MCZ 

M320964 

MCZ 

M320965 

MCZ 

M320967 

MCZ 

M320968 

MCZ 

M320971 

MCZ 

M382341 

MCZ 

M382364 

MCZ 

M382380 

MCZ 

N141771 

NMNH 

N 152408 

NMNH 

N592677 

NMNH 

N133176 

NMNH 

N150528 

NMNH 

N152520 

NMNH 

N152527 

NMNH 

N153056 

NMNH 

N153057 

NMNH 

N181808 

NMNH 

N 184442 

NMNH 

N184833 

NMNH 
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